Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  814 1467 (N31E16 L= 38)   437     85        1 ( 770/1440- 853/1484)   9.5 E19-E13 N29-N32
      2 1829 1330 (N22W67 L=121)   853    131        0 (1771/1284-1871/1368) -26.6 W60-W73 N18-N24
      3 1511 1311 (N21W33 L= 88)   266    101        0 (1469/1273-1541/1338) -38.7 W30-W35 N18-N23
      4  388 1219 (N13E45 L=  9)   858    251        3 ( 352/1113- 442/1317) -55.3 E49-E39 N06-N19
      5 1096 1252 (N17W04 L= 58)   404    101        0 (1044/1233-1136/1271)   9.5 W00-W06 N16-N18
      6 1764 1220 (N15W55 L=109)   333     63        1 (1732/1202-1793/1244) -18.4 W51-W57 N14-N16
      7  648  942 (S03E24 L= 31)   259     84        0 ( 638/ 899- 659/ 986) -90.0 E24-E22 S05-N00
      8 1736  846 (S08W51 L=105)   303     65        0 (1710/ 821-1760/ 876)  45.0 W48-W52 S09-S05
      9  221  828 (S12E61 L=353)   991     92        5 ( 183/ 790- 249/ 857) -29.7 E66-E57 S14-S09
     10  736  604 (S24E19 L= 35)   284     93        1 ( 700/ 572- 770/ 638) -38.7 E21-E16 S26-S21
     11 1401  430 (S36W31 L= 86)   851    181        3 (1373/ 350-1424/ 518) -81.3 W27-W36 S42-S29
  

    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