Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1632 1625 (N28W67 L=230)   494     83        1 (1605/1605-1670/1652) -53.1 W64-W69 N25-N29
      2 1397 1470 (N25W35 L=198)  1591    322        1 (1257/1373-1481/1580) -71.6 W30-W39 N17-N35
      3 1266 1190 (N12W18 L=181)   467    207        0 (1207/1133-1355/1217)   0.0 W13-W23 N09-N13
      4 1617  920 (S12W36 L=198)   415    173        0 (1589/ 887-1653/ 970) -90.0 W33-W38 S13-S08
      5 1381  935 (S05W20 L=182)   413    135        1 (1327/ 913-1454/ 957) -32.0 W16-W24 S08-S03
      6 1652  748 (S23W37 L=200)  2147    532        2 (1499/ 584-1788/ 956)  30.7 W23-W51 S30-S13
      7 1493  844 (S13W26 L=188)   359    103        0 (1469/ 796-1510/ 894)  59.0 W23-W27 S15-S10
      8  439  864 (N07E41 L=122)   248     40        1 ( 418/ 845- 460/ 882)  14.0 E42-E38 N05-N07
      9 1740  695 (S30W46 L=209)  1296    246        1 (1674/ 606-1815/ 788)  31.0 W39-W54 S34-S25
     10  620  571 (S14E35 L=128)   311     73        1 ( 585/ 547- 654/ 586)   0.0 E37-E32 S15-S12
     11  549  467 (S20E45 L=118)   358     75        1 ( 531/ 429- 570/ 492)  29.7 E48-E41 S22-S17
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y